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A new definition of a dipole glass parameter for ferro-antiferroelectric solid mixtures is suggested. It is con- 
structed on the nearest neighbours pair correlation functions for interacting dipole momenta. The behaviour of 
the dipole glass parameter is calculated and discussed. 
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1. Introduction 

The problem of a dipole glass state in frustrated ferroelectrics is difficult both in experimental detec- 
tion and for theoretical interpretation. It still remains actual, because this state of non-ordering system 
at low temperatures is widely observed in different materials, for example in dielectrics with non-central 
ions, relaxors, mixed ferro-antiferroelectric compounds . Dipole glass state is a characteristic one 
for structurally complex systems and corresponds to non-equilibrium thermodynamic phases in them. 
The most important property of the system and a simple criterion for detection of the process of dipole 
glass transition is a dispersion of dielectric characteristics. Therefore, in experimental investigations, the 
dynamical methods for the study of real and imaginary parts of dielectric susceptibility behaviour under 
high and low frequencies are used. Naturally, this state does not depend on the external field effect and 
must be observed in a static case. 

On the other hand, due to a competition between different types of interparticle interactions in a com- 
plex crystalline matrix there is a possibility of different local states coexisting with their non-equilibrium 
occupation. Those states correspond to a local minima of the free energy of the system with high poten- 
tial barriers between them. In the thermodynamic limit (N — ► oo, V — ► oo, Nl V = const), the relaxation 
time for transition from some state to another one tends to infinity. As far as the system is close to such a 
local minimum, it remains there for a long time. These systems are not ergodic and their theoretical de- 
scription is quite complicated. The dispersion of dielectric characteristics of the systems in a dipole glass 
state is caused by the stimulated transitions from a certain energy state to another one. So, the dynamic 
behaviour of the investigated system can be illustrated by their microscopic structure, or in other words, 
by their static correlation functions. 

The experimental study of dipole glass phases in frustrated ferroelectrics is presented mostly 
in papers devoted to mixed ferro-antiferroelectric compounds based on the KDP-type crystals 
[K„(NH 4 )i-„H 2 P04, Rb„(NH 4 )i-„H 2 P04, Rb n (NH 4 )i_„D 2 P0 4 , Rb„(NH 4 )i-„H 2 As0 4 ]. The most fun- 
damental of them, in our opinion, are dielectric investigations Ifjjllll. In all of them, the temperature 
behaviour of real (e'j and imaginary (e") parts of dielectric permittivity at low and high frequencies looks 
similar for concentrations 0.2 < n < 0.8. Namely, the existence of smooth (e'j or sharp (c") peaks at tem- 
peratures about 50 K indicate some ordering processes in the system. Since no long-range ordered phases 
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in this range of concentrations were observed, only an abnormal increase of short-range correlations is 
possible there. This state was called a dipole glass phase. 

The EPR NMR (3 and X-ray (Hi structural investigations have shown the existence of signif- 

icant local fields in ferro-antiferroelectric mixtures at low temperatures. Those fields are connected with 
strong short-range correlations between charged particles (protons). 

The theoretical study of mixed ferro-antiferroelectric systems is based on the use of a replica method 
for calculation of the basic thermodynamic functions. This method was first suggested in the theory 
of diluted magnetic or mixed ferro-antiferromagnetic systems 1 16-20]. As for mixed hydrogen bonded 
ferro-antiferroelectric compounds, the Ising model with a random transverse and longitudinal fields 
was explored. The Gaussian distribution for random interactions and the corresponding fields was used 
in 1 21,|22|]. A cluster theory of mixed KDP-type ferroelectrics, including the effect of piezoelectric coupling, 
was proposed in (23^250. The phase diagrams, the behaviour of thermodynamic functions, long-range and 
dipole glass order parameters were obtained and discussed. Most of those results are in good agreement 
with the experiment. 

In our previous papers litj - Eill . a new approach to the study of hydrogen bonded ferro-antiferro- 
electric solid mixtures of Rb„(NH 4 )i_„H2As04 (RADA) type was suggested. It is based on the concept of a 
decisive role of microscopic structure of a system (single and pair correlation functions) in the formation 
of a dipole glass phase. The cluster mean field approximation with replica method was used for their 
calculation. So, the "concentration-phase transition temperature" phase diagram was obtained in good 
agreement with the experimental one. The areas of the existence of dipole glass phases were interpreted 
as multiplication of pair correlation functions at a certain concentration for one of the components at 
relatively low temperature. Since the pair correlation functions characterize a short-range order in the 
system, the dipole glass state is determined by different possibilities of short-range arrangement of inter- 
acting particles, i.e., by their stable and metastable states all over the sample. 

As for a strict definition of a dipole glass parameter, there appear to be some problems. Usually 1 16] it 
is defined as a mean value of the dipole momentum square. However, in this case, the multi-stage charac- 
ter of dipole glass phase and the real interparticle correlations are not taken into account. Moreover, the 
temperature behaviour of dipole glass parameter does not react adequately upon passing the "freezing" 
point (the point of the dipole glass phase appearance). This parameter remains non-equal to zero even 
in the paraelectric phase where no glass state can be observed. A more appropriate definition of this 
parameter and the static interpretation of a dipole glass phase will be presented in the present paper. 

2. Model. Analytical results 

It is well known that 1 29] the ferroelectric RbH2As0 4 and antiferroelectric NH4H2ASO4 single crystals 
are isomorphic with tetragonal lattice space group symmetry in paraelectric phase. Below 110 K 
in the first crystal and below 213 K in the second one, a ferroelectric and antiferroelectric ordering are 
correspondingly observed. The axis of the orientation of dipole moments in RbH2As04 is directed along 
OZ, but in NH4H2ASO4 it is directed along OX (07). 

Solid mixtures of both crystals due to frustration of ferro- and antiferroelectric types of interactions 
demonstrate a possibility of a new "ordered" dipole glass phase appearance. In this phase, no type of 
long-range ordering is preferential and only short-range correlations between particles remain. Dipole 
glass phase is realized for intermediate concentration (0.50 < n < 0.80) and at relatively low tempera- 
tures (below 100 K). The physical reason for the existence of these phases are numerous arrangements of 
nearest neighbours (particles) in a disordered ferro-antiferroelectric system. 

We will regard a regular KDP-type crystalline lattice, each site of which is randomly occupied by 
rubidium atoms or ammonium groups, that is by Rb or NH4. The effective Hamiltonian of this mixed 
system is taken in the form 1 28, 30]: 

H = -\L E \^ m 'n i n j S z im S) m , + UT j m \l-n i ){\-n j )S im S) m \ 

i,j mm' 

+ ^EE [£cos0%Sf m + Esin0(l - hi)S x im \ , (2.1) 

£ \ m 
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where both V™ m and U™ m are the intensities of interaction between m-th dipole particle in i-th site 
with m'-th dipole particle in 7-th site when dipole moments are oriented along z-axis and along x-axis, 
correspondingly. S z and S x are projections of the unit classical vector S, the site occupation operators hi 
have the following eigenvalues: 

1, when i-th site is occupied by Rb atom, ^ 
0, when i-th site is occupied by NH4 group. 

Here, E is the external electric field, 8 is the angle between z-axis and the direction of this field. 

To take into account the random distribution of dipoles oriented along z-axis (Rb-AsC>4 groups) and 
dipoles oriented along x-axis (NH4-ASO4 groups) we have used the replica method 1 18] for the configu- 
ration averaging of thermodynamic functions. Since there are two types of dipoles having ferroelectric 
and antiferroelectric orientation, we must correspondingly introduce two types of replica variables o k im 
for S z im and fjv for S x im . The following property of a binomial random variable %c will be used. If the 
probability of Xc is Pr{x c - r } - (r) P r q cr > then the expectation value of exp{ax c } is equal to [pe a + q) c , 
which behaves like 1 + cln(pe a + q), when the analytic continuation for c — » is performed. In the con- 
text of the investigated system, the p and q are as follows: p — n, q = 1 - n, n is a concentration (density) 
of Rb atoms. 

Since the unit cell for KDP-type crystal contains two formula units (there are also two sublattices 
below the phase transition point 1 29]), the numbers m, m' run from 1 to 2. So, we consider the next form 
for the interactions between all particles in the crystalline lattice: the strict two-particle interaction for 
particles belonging to the same crystalline lattice site and the self-consistent mean field approximation 
for particles from different sites, the so-called two-particle cluster approximation filll . 

With the accuracy up to the second order of symmetric replica expansion and using two-particle 
cluster approximation, the following expression for the free energy of the system investigated has been 
obtained: 

F = -iln{4[e" Al +e Al cosh(2£i+/3£cose)] 

x [e A2 cosh(/3£sine) + e" A2 cosh(2£ 2 )]} + -C. (2.3) 



Here, 



Ai = h+h{v 2 f -V 2 af+ \\ + {Vf-V af f\g z -Ug x ), 

A 2 = (tanh pV-h)U + J 2 {u}-Ul f +\u 2 + {Uf-U af ) 2 ]g x -Ug z }, 

By = \hV f + A] 2 [vjP 2 -W f U af q 2 ^\P, 

B 2 = [[tznhpV-h)U af +4J 2 (ul f q 2 -3V f U af P 2 )]q, (2.4) 

C = /iV / P 2 +(tanhj3V-/i)f> B/ <f 2 +Y \-l - U + {g z - g x f 

+ (?/ + %) 2 + ~ Vaff g 2 z + {Of + U af f + {U f + U af ) 2 g 2 x + 12 [V f P 2 - U af q 2 f] , 



v - Vf v Vaf fi Uf fi Uaf 
v f-y v «f-—> u f~y u *f- — ' 



N , , 1 JV 



V = V} 2 , U = U} 2 . (2.5) 
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oo 

/ 1 = X(-l) t 2ln(l-n + ne- 2 ' ,vt ) l 

oo 

J 2 = £(-l) f 4fln [l-n+ne~ 2fivt , (2.6) 
t=i 

/3 = (fcT) -1 , fc is the Boltzmann constant, T is the absolute temperature. 

The form of ]\{n, T) and Jzin, T) functions is determined by the binomial random distribution of 
different types of "particles" (Rb or NH4) and the replica method of configuration averaging. ]\{n, T) is 
proportional to the concentration n, but ]%{_n, T) reflects the correlation between both subsystems and 
naturally tends to zero in the limits n — ► and 

Ferroelectric P, antiferroelectric q order parameters and g z , g x nearest neighbours pair correlation 
functions are defined in the following way: 

P = < on >, q = <Sn>, g z = <<JiiOi 2 >, g x = <$n$n>, (2.7) 

where < . . . > denotes thermal averaging for a given distribution of z- and x-oriented dipole momenta of 
Rb-As0 4 and NH4-ASO4 groups, respectively, whereas a bar denotes a stochastic averaging over different 
random distributions. 

From the principle of stationarity for the free energy H2.3t under order parameters P, q and both pair 
correlation functions g z and g x (they determine a structure of nearest neighbours particles) the set of 
equations is obtained: 

sinh(2Bi +/3£cos0) 
cosh(2Bi + /3£cos0) + e" 2A i ' 

sinh(2B 2 ) 
cosh(2B 2 ) + e 2Az cosh(/3£sin0) ' 
-1 + e 2Al cosh(25i + /3£cos0) 
1 + e 2A i cosh(2Bi + fiE cos 0) ' 
cosh(/3£sin0) - e~ 2A2 cosh(2£ 2 ) 

cosh(/3£sin0) + e~ 2A * cosh(2B 2 ) ' (2 ' 8) 

Spontaneous values (in zero external field E) of P, q, g z ,gx parameters are limited by + 1. So, naturally, 
they may be interpreted as a set of order parameters of the system under investigation. 



P = 

q = 

gz = 

gx = 



3. Phase diagram. Dipole glass parameter 

Computing the equations <2.8> one can find the dependencies of P, q, g z , g x parameters on tem- 
perature T and concentration n. Those dependencies form a phase diagram of the investigated ferro- 
antiferroelectric mixed system in the plane n - T. The principal outline of a phase diagram for RADA 
system was obtained in the paper 1 28], but due to its considerable role in the analysis of the mixed system 
behaviour, we reproduced it here (figure[l). 

The main features of the phase diagram are: 1) the existence of two non-crossing areas for ferro- 
electric (F) and antiferroelectric (AF) phases; 2) the existence of pure dipole glass phase at intermediate 
concentration and relatively small temperature (this phase is divided into several (G3-G9) subphases ac- 
cording to the order of the chaotic state); 3) the possibility of mixed ferroelectric-dipole glass (F-G) and 
antiferroelectric-dipole glass (AF-G) phases existence. It should be noted that the upper dipole glass phase 
G3 is rather a precursor of the real dipole glass phase appearance. This conclusion is in good agreement 
with the experiment (HIH. 

The difference between RADA phase diagram in 1 28] and the one presented here is the line divid- 
ing the paraelectric phase and the antiferroelectric phase (P-AF). In figure [T] the upper line represents 
the points of stable non-zero root arising for antiferroelectric order parameter q while in 1 28] the cor- 
responding line indicates the points of its bifurcation (lower line in figurejl). It should be noted that the 
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only two fitting points were used to build figure [T] The first one is the antiferroelectric phase transition 
point temperature for pure NH4H2ASO4 crystal, and the second one is the ferroelectric phase transition 
point temperature for pure RbH2As04 crystal. 

The dipole glass phase is interpreted as a possibility of different arrangements for the particles of 
nearest neighbours. This fact is reflected by numerous solutions for pair correlation functions g z and g x 
at different values of temperature and concentration. When the temperature decreases, the dipole glass 
state becomes more complicated. 

In this context, the behaviour of correlation functions g z , g x is characteristic of the dipole glass struc- 
ture. However, there are a some peculiarities connected with their relation to dipole glass state parame- 
ter. Usually 1 34], this parameter is introduced to ferroelectrics theory by analogy with spin glass param- 
eter in ferromagnetic systems 116], named Edwards and Anderson parameter. Edward and Anderson 
parameter is defined as a second moment of magnetic order parameter. This moment is maximum for 
low temperature and vanishes at relatively large temperature T g named glassy temperature. As for ferro- 
electric systems, the situation is not that simple. The recent investigations 1 25] of mixed hydrogen bonded 
ferro-antiferroelectric systems have shown some restriction to using the Edwards and Anderson param- 
eter for a complete definition of a dipole glass state. The second moment of ferroelectric order parameter 
does not vanish even for very high temperatures when no glass state can exist. 

Pair correlation functions g z and g x , which are analogous to the second moment of order parameter, 
are nonzero even in paraelectric phase P (see figure [I) and tend to zero for very high temperatures 
(theoretically limj-— oogz = lim^-. 00 gx = 0). According to our concept, the dipole glass phase is a state 
with numerous types of short-range correlations between particles, or numerous numbers of solutions 
for a set of equations <2.8> . The intensity of correlations (values of g z and g x ) is important but not decisive 
for dipole glass phase characterization. Dipole glass phase arises only at the moment when the number 
of different solutions for g z or g x changes from 1 to a bigger number. Such a situation is observed on the 
upper line of G3 area (figure[l}. As we have calculated (the phase diagram has been built on this basis), 
there is only one solution for (g z , g x ) pair in the paraelectric phase (P), three solutions for {g z , g x ) pair 
in the G3 area, five solutions for (g z , g x ) pair in the G5 area and so on. With an increase of the number of 
different solutions for g z , g x , a bigger number of different types of correlations of the nearest particles 
are observed. As a result, a more complicated arrangement of dipoles all over the crystal takes place. This 
situation is accompanied by a more chaotic distribution of correlations and by the formation of a dipole 
glass phase. When temperature decreases, a higher stage of dipole glass phase is realized. Thus, dipole 
glass parameter should take into account two factors: the number of solutions for g z , g x parameters 
and their absolute values. It is obvious that the unique (trivial) solution for g z , g x , which remains in the 
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paraelectric phase must be eliminated. 

Therefore, since both g z and g x functions form a dipole glass parameter G, we propose the next 
expression for its calculation: 

\ m-l , 



G = 



2M 



r 1 

£ [n\gM\ + a-n)\g?\\. 



(3.1) 



where , g^' are the roots of a set of equations (2.8) , m is the number of these roots for any temperature 
and concentration, No is the maximum number of roots, which are taken into account (in the presented 
calculations Mq = 9), n is the concentration. 

The "partial" dipole glass parameters G z and G x : 



j m—l 

1 m—l 
™0 (=1 



(3.2) 



describe the power of dipole glass state intensity according to z - z and x-x correlation of Rb-As0 4 and 
NH4-ASO4 dipoles, correspondingly. 
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Figure 3. Temperature dependencies of the dipole glass parameters G for different concentrations. 
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The temperature dependencies for G z and G x at different concentrations are presented in figure [2] 
and for G they are presented in figure [3] The most characteristic feature in the behaviour of G z , G x and 
G are their step-like dependence. The multigraded character of a dipole glass parameter is caused by the 
complex nature of dipole glass phase to exactly change the number of roots for g z , g x correlation func- 
tions (see figure [l}- The position of "jumps" for G z , G x , G is determined by the lines dividing the areas 
of P-G3, G3-G5, ... for "pure" dipole glass phases and F-G3, G3-G4, . . . ; AF-G3, G3-G5, ... for the mixed or- 
dered and dipole glass phases. The existence of those lines is determined by the finite number of different 
roots of equations d2.8> in every area of temperature and concentration. Theoretically this number is in- 
finite (32|,[35|], so the step-like behaviour of dipole glass parameter in these limits may transform into a 
quasi-continuous curve. 

4. Conclusions 

A new approach to the calculation of the dipole glass parameter in mixed ferro-antiferroelectric sys- 
tems has been developed. It is based on the conception of a decisive role of the pair correlations of 
the nearest neighbours particles in a crystalline lattice for the dipole glass phase formation. The cor- 
responding correlation function has been calculated in the second order of replica expansion and using 
the cluster approximation. This microscopic approach leads to the basic physical statement: the dipole 
glass phase is a polystate arrangement of the nearest neighbours particles with the possibility of non- 
equilibrium occupation of them. In such a way, the proposed definition differs from the well known 
Edwards and Anderson spin glass parameter in the quenched ferromagnetic system theory. 

The calculated dipole glass parameter takes into account not only the intensity of interparticle correla- 
tions but also the multiplicity of a different type of those correlations. The latter corresponds to different 
possible types of short-range order in the mixed system and usually describes a metastable states. The to- 
tal dipole glass parameter is formed both by ferro-ferro and antiferro-antiferro types of correlations. The 
step-like behaviour of the dipole glass parameter reflects the non-uniform structure of the dipole glass 
phase and the dependence of the intensity of the disorder on temperature and concentration. When the 
number of different roots for pair correlation functions increases, the temperature dependence of this 
parameter becomes more smooth. 

The obtained phase diagram "temperature-concentration" demonstrates the possibility of the exis- 
tence of both pure ferroelectric, antiferroelectric, dipole glass phases and mixed ferroelectric-dipole 
glass, antiferroelectric-dipole glass phases. This fact was recently confirmed by numerous experimental 
investigations. The precursors of the dipole glass phase which appears at rather high temperatures are 
interpreted as the state with the finite number of different types of interparticle correlations, in counter- 
balance to paraelectric state with only one type of those correlations. 
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noBefliHKa napaiueTpa fli/mo/ibHoro cxsia pji* 
cerHeTO-aHTi/icemeToe/ieKTpi/iHHi/ix TBepfli/ix cyMiweft 

M.A.Kopi/iHeBCbKi/ibPE, B.E.Co/iOB'flhP 

iHCTHTyr c|)i3HKW KOHfleHCOBaHnx cucreM HAH YKpaTHM, By/i. I. CBCHL^ii^bKoro, 1, 7901 1 /lbBiB, YKpaiHa 
Hai^iOHa^bHwii yHiBepcuTeT "ZlbBiBCbKa no/iiTexHiKa", syn. C. EaHflepw 12, 79013 /lbBiB, YKpaiHa 
iHCTMTyT 4>i3HKW, LUeijiHCbKHti yHiBepcmeT, syn. Be/ibKono/ibCKa 15, 70451 LUei^iH, llo/ibma 

3anponoHOBaHO HOBe 03HaHeHHa napaMeipa flwno/ibHoro cioia fl/ia cerHeio-aHTucerHeToe^eKTpuHHwx TBep- 
py\x cyMiiuePi. L(ePi napaMeTp 6yflyeTbca Ha napHwx Kope/ifmmHwx tfiyHKi4inx B3aeMOfliK)4wx flwno/ibHwx mo- 
MeHTiB nacTHHOK Haw6/inx4nx cycifliB. Po3paxoBaHO Ta floaiiflxeHO noBefliHKy napaMeipa fluno/ibHoro cua, 
a TaKOx c|)i3W4Hi B/iacTHBOCTi (f)a3n flwno/ibHoro CK/ia. 

K/i h) h o b i c/ioBa: cerHeToeneKTpi/iKM, amiAcerHeioeneKipiAKvi, 3Miiuat-ii CHCTeMM, Kope/ifujiiiHi cpyHKLfif, 



13703-8 



